Magnetoinductance of Josephson junction array 
with frozen vortex diffusion 
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The dependence of sheet impedance of a Josephson junction array on the applied magnetic field 
is investigated in the regime when vortex diffusion between array's plaquettes is effectively frozen 
due to low enough temperature. The field dependent correction to sheet inductance is found to be 
proportional to /ln(l//), where / <C 1 is the magnitude of the field expressed in terms of flux 
quanta per plaquette. 
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I. INTRODUCTION 

A two-coil mutual-inductance technique [1,2] has been 
introduced by Fiory and Hebard [1] for measurement of 
(linear) frequency dependent complex sheet impedance 
Za(u>) [or, equivalently, sheet conductance Gn(u>) — 
1 / Z a (w)] of thin superconducting films [3] . This approach 
has also proved itself very useful for investigation of ar- 
rays of weakly coupled superconducting islands (Joseph- 
son junction arrays) in external magnetic field (for a re- 
view, see Ref. [4]). In such systems the nature of the 
response is strongly dependent on the magnitude of the 
applied dc field, which is convenient to discuss in terms 
of the ratio / = <&/$() of the magnetic flux <i> penetrating 
each array plaquette to the superconducting flux quan- 
tum $o- 

For an integer (for example, zero) / a Josephson junc- 
tion array can be described by the two-dimensional XY 
model and with decrease of temperature experiences the 
Berezinskii-Kosterlitz-Thouless transition [5,6] into a su- 
perconducting phase. In this phase all vortices are bound 
in neutral pairs and the sheet impedance in the limit of 
low frequency becomes purely inductive: Z a (oj) ~ iujL a , 
the effective sheet inductance L a being inversely propor- 
tional to the superfluid density. A deviation of / from an 
integer value (by 8f) introduces a finite concentration of 
unbound vortices c (in dimcnsionless units, that is per 
lattice plaquette, c = \Sf\), as a consequence of which 
the response is strongly changed and in the low frequency 
limit becomes purely dissipative: Zu(lS) w cRy (the fre- 
quency independent constant Ry being associated with 
a contribution from a single vortex). 

Analogous behavior can be expected to take place in 
the vicinity of every rational /, if the temperature is 
lower than the (discontinuously dependent on /) tem- 
perature of the phase transition associated with freezing 
of the field-induced vortices into commensurate pattern 
[7]. As a consequence, a measurement of the real and 
of the imaginary components of Z a (w) as a function of / 
(at fixed to and reasonably chosen temperature) produces 
strongly oscillating curves [8-12], the well developed dips 
on which correspond to the superconducting states. 



With further decrease of temperature the dissipative 
component of Z a (w) becomes very small and the oscilla- 
tions of the inductive component much less pronounced 
[8,9,12]. This is related with suppression of vortex diffu- 
sion at temperatures for which the rate of a thermally 
activated tunnelling of a vortex between neighboring lat- 
tice plaquettes (which requires to overcome a well defined 
barrier [13]) becomes much smaller than the frequency of 
the measurement. 

The present work is devoted to investigation of a 
Josephson junction array response in this particular 
regime, when the relation between the frequency and 
the temperature allows to consider all vortices as frozen 
in lattice cells which they occupy. In terms of a simple 
model [4,13] treating a vortex in a proximity junction 
array as an overdamped point particle moving in the ex- 
ternal potential (imposed by the structure of the array) 
this corresponds to the case when each of the vortices 
is confined to oscillate within the limits of a particular 
minimum of the potential, which then can be replaced by 
a harmonic one. In the framework of such description the 
correction to sheet inductance (per vortex) would have 
a finite value, inversely proportional to the curvature of 
the potential. 

With the help of a more straightforward calculation 
(based on the reduction to the equivalent electric circuit) 
we show that the single vortex contribution contains log- 
arithmic divergence and, therefore, for small magnetic 
fields (that is small vortex concentration, /< 1) the cor- 
rection to Lp 1 is proportional to /ln(l//). The results 
can be of interest in relation with sheet impedance mea- 
surements for small / at intermediate frequencies [9,10]. 



II. THE HAMILTONIAN AND THE VORTICES 



In absence of external magnetic field a regular Joseph- 
son junction array can be described by the Hamiltonian 



(1) 



K) 



1 



where p n is the phase of the order parameter of n-th 
superconducting island. In the case of a square lattice 
the variable n can be chosen in the form of the vector 
n = (n x ,n y ) with integer components n x and n y . The 
summation in Eq. (1) is performed over all pairs (nn') of 
coupled islands and 

j = A /c (2) 

is the Josephson coupling constant, which is assumed to 
be the same for all junctions, I c being the critical current 
of a single junction. The form of Eq. (1) assumes that the 
coupling is weak enough, so the magnetic field created by 
the currents flowing in the array can be neglected. 

Variation of the Hamiltonian (1) with respect to ip a 
shows that the minimums of H are achieved when the 
variables {p> n } satisfy the current conservation equations 

£ J nn'=0, (3) 
n' 

where 

Inn' = Ic SUl(yV - (p n ) (4) 

is the current from the n-th to the n'-th island, defined 
only for the pairs of the islands connected by the junction. 
The simplest solution of the Eqs. (3)-(4) is the trivial so- 
lution 



Continuous approximation can be also used for esti- 
mating the energy of a vortex 

E v = J5>-cos(rii (9) 

nn') 

because the integral 

Ev ~~2-f j dxd V&<P v ) 2 , (10) 

to which the lattice sum of Eq. (9) is reduced in the 
framework of this approximation, diverges at the upper 
limit: 

E v w 7r 7 Jln7V . (11) 

Here N is the linear size of the array and 7 is the numer- 
ical constant depending on the structure of the lattice 
(for square lattice 7=1 and for triangular one 7 = \/3). 
Accordingly, the main contribution to Eq. (9) is coming 
from the large scales, where it can safely be replaced by 
its continuous form (10). As a consequence of this diver- 
gence, at low temperatures all vortices, which appear as 
thermal fluctuations, are bound in neutral pairs. On the 
other hand, a finite concentration of vortices of the same 
sign can be induced by application of external magnetic 
field in perpendicular to array. 



p n = const (5) 

corresponding to the global minimum of H and the ab- 
sence of any currents. 

A vortex is a local minimum of H, in which on going 
along any closed loop surrounding the vortex core (which 
can be associated with a particular plaquette of the lat- 
tice) (p n changes by 2irs, where s = ±1 is the topologi- 
cal charge of a vortex. The form of this solution implies 
the presence of persistent currents circulating around the 
vortex core. 

Away from the core <p a changes slowly and Eqs. (3)-(4) 
can be linearized to give: 

- ¥>n) = , (6) 

n' 

where the summation [like in Eq. (3)] is performed over 
the nearest neighbors of n. In continuous approximation 
Eq. (6) is reduced to 

VV = , (7) 

which allows to conclude that at large distances from the 
vortex core the spatial distribution of p> is given by 

(fi v (x, y) « s arctan — + const , (8) 

x 

where x and y (r 2 = x 2 + y 2 3> 1) are the coordinates 
(in lattice units) counted from the core. 



III. THE EQUIVALENT ELECTRIC NETWORK 

Expansion of the Hamiltonian (1) up to the second 
order in deviations 

<fyn = Pn - (12) 

of the variables tp n from their values <p„ ^ in some of the 
minima of H gives 

H<® {Sip} = H{tpM + Sp} - H{p^} 

~^Jnn'(<^n-<W) 2 , (13) 
nn' 

where 

Jnn' = JCOS^-^]. (14) 

Since the deviation of the current SI nn > from its value in 
the extremal solution 7 nn / = Iq sm[pn ^ — <p$] m l mear 
approximation is given by 

Sinn' - Ic COS^ - p^)}(Sp n - Sp n ,) , (15) 

Eq. (13) can be rewritten as 

H{2) = \ E inn'O^nn') 2 > ( 16 ) 
(nn') 
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where 



and 



inn' 



u = 



L 



r (o) (0)1 
cos[<^„ ; - ip n /\ 



1 2 



J 2el r 



(17) 



(18) 



The form of the Hamiltonian (16) allows to conclude 
[14,15] that in the harmonic approximation the array be- 
haves itself with respect to additional (external) current 
as the network formed by the inductances L nn > = L n > n 
defined by Eqs. (17). At finite frequencies it is also neces- 
sary to take into account that each of these inductances 
is shunted by the resistance R nn ', so the complex con- 
ductance G nn < (w) of each network link is given by 



G„„'M 



1 



1 



(19) 



In the case of a proximity coupled array, which is ex- 
plicitly considered below, the shunting resistance R nn ' 
is determined mainly by the conductivity of the under- 
lying metallic substrate, and, therefore, can be assumed 
to be frequency independent and the same for all the 
links: i? n n' = Ro - In absence of vortices tp a = const, and, 
therefore, the conductance of all the links is the same: 



Gnn'M = G (oj) 



1 1 

iuLq Rq 



(20) 



The distribution of the currents I nn '(w) in such net- 
work has to be found by solving the current conservation 
equations of the form (3) with 



Inn>{u)=G nn ,(w)(y n -V n ,) , 



(21) 



where V n is the amplitude of the time dependent electric 
potential 



V a (t) = V n exp(iut) 



(22) 



on the n-th superconducting island. Application of the 
time dependent potential difference 



V(t) = Vexp(iwt) 



(23) 



for example in the x direction to the square network 
N x N formed by equivalent elements, leads to the dis- 
tribution of the potentials 



const 



and the currents: 



G {lo)V/N forn' = n + e x 
for n' = n + e y , 



(24) 



(25) 



[where e x = (1,0) and e y = (0, 1)] corresponding to the 
total current in the chosen direction given by 



J(w) = G (tu)V . 



(26) 



This means that the sheet conductance G D (ui) of a uni- 
form square network coincides with the conductance of a 
single link: 



G n (w) = G (w) . 



(27) 



In case of a uniform triangular network Gn (oj) = 
V / 3Go(w) and, like in the case of a square network, the 
sheet conductance does not depend on the direction. 

In the following it will be convenient to decompose 
Gn (w) into the real and the imaginary parts as 



G n (u;) 



1 



iu>Ln(u>) i?rj(w) ' 



1 



(28) 



where the effective sheet inductance L a {oj) and the ef- 
fective shunting resistance Ru (u>) are the real functions 
of uj. In particular, Eq. (27) corresponds to L n (u>) = Lq 
and i?o(w) = Rq- 



IV. THE SINGLE VORTEX CORRECTION TO 
CONDUCTANCE 

In order to find the correction to frequency dependent 
sheet conductance related with the presence of a vor- 
tex (which is assumed to be frozen in a particular array 
plaquette) one has to consider a network in which the 
distribution of ip n and, therefore, of G nn < is assumed to 
correspond to vortex configuration. For non-uniform ip n 
it is convenient to rewrite the expression for the total 
current (in the x direction) in a square network as 



N 
1 



E G 



V 



(29) 



where 



Pn = G n>n + e — G 



ILU 



(30) 



and v n parametrizes the deviation of V n from its value in 
a uniform network: 



t/ n x + v n 

V n = — — V + const 



(31) 



and has to be found by solving the current conservation 
equations, which for V n of the form (31) can be rewritten 
as 



^G nn '(o>)(u n - V) = Pn 



(32) 
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Comparison of Eq. (29) with Eq. (26) shows that the 
correction to Go (w) induced by a non-uniform distribu- 
tion of <p n can be split into two contributions, the first of 
which has the form of a frequency independent correction 



(^□ 1 )l--^E( i 1 -in, 1 n + e, 



N 2 2J 



(33) 



and turns out to be proportional to the energy E 
(counted from the ground state energy) of the considered 
configuration. In the case of a single vortex this energy 
is given by Eq. (11) and, accordingly, the expression for 
(ALq 1 )! contains the logarithmic divergence. 

The form of this correction corresponds to formal av- 
eraging of the superfluid density. In conjunction with the 
absence of a frequency dependence this allows to conclude 
that (ALq 1 )! has nothing to do with vortex oscillations 
in a potential minimum and should be associated with 
suppression of superfluid density. 

The second correction 



(AG D ) 



1 



(34) 



is characterized by a more complex frequency dependence 
and reduces to correction to L^ 1 only in the limit of 
u> — > 0. It is not hard to show that in the case of a sin- 
gle vortex the lattice sum in Eq. (34) [in contrast to the 
lattice sum in Eq. (33)] is not divergent at large scales. 

The behavior of P n away from the vortex core can be 
found from the continuous approximation, in the frame- 
work of which 



d 



iuL dx 



cos 



dip 
dx 



1 



d 



'dip 

2iujL dx V dx 



For (p(x, y) of the form (8) this gives 



iuiLo (x 2 + y 2 ) 3 



£ -(35) 



(36) 



Therefore, if v n decays with the increase of the distance 
from the vortex core (as naturally one can expect it 
to do), the lattice sum in the right-hand side of Eq. 
(34) will be convergent and the behavior of correction 
to ^□ 1 ((jj = 0) will be dominated by the divergence of 
the contribution (AL^ 1 )! discussed above. 

The behavior of v n away from the vortex core can 
be found by replacing in the left hand side of Eq. (32) 
G nn '(w) by Go(w) = l/iujL n + 1/Ro, after which it is 
reduced to 



— (1 + iLuL /R )V 2 v = iLuL P(x, y) . 



(37) 



The solution of Eq. (37) allows to find that away from 
the core v{x,y) decays as 



v(x,y) 



1 



iujL /R 
x 



(38) 



4(x 2 + y 2 ) 



ln(*'+ir) 



2U/2 



x(3y 2 - x 2 



16(x 2 + y 2 ) 2 



which confirms the convergence of the lattice sum in Eq. 
(34). 

The frequency dependence of (AG n )2, related with the 
frequency dependent factor in Eq. (38) is consistent with 
what one expects from the driven oscillations of an over- 
damped particle in effective harmonic potential. Quite 
naturally (AL^ 1 ^ = lim^^o iw(AGn)2(w) [as well as 
(ALg 1 )!] does not depend on the value of shunting re- 
sistance Rq. 



V. RESULTS AND DISCUSSION 

Thus we have found that at frequencies for which the 
diffusion of vortices is effectively frozen the main cor- 
rection to the sheet inductance of a Josephson junction 
array comes from the suppression of the superfluid den- 
sity rather then from the displacement of vortices in the 
effective potential and contains logarithmic divergence. 

In the presence of a small concentration of more or less 
uniformly distributed vortices induced by the application 
of weak magnetic field (with / <C 1) the divergence in the 
expression for the vortex energy is cut off at r ~ / _1 ^ 2 
(instead of at r ~ N), and the correction to L^ 1 can be 
rewritten as 



4 / 



(39) 



The same expression is also valid for other periodic lat- 
tices, for example triangular. When the screening effects 
related with the self-induced magnetic fields of the cur- 
rents in the array are taken into account, the logarithmic 
factor in Eq. (39) has to saturate (with decrease of /) 
when the typical distance between vortices becomes of 
the order of the magnetic penetration depth [16,17] of 
the array. 

Earlier the magnetoinductance of a Josephson junction 
array in the regime of frozen vortex diffusion has been 
investigated for the case of a fractal array with the struc- 
ture of the Sierpinski gasket [15,18]. With the help of the 
recursive calculation using the self-similarity of the Sier- 
pinski gasket it has been found [15] that the field depen- 
dence of the correction to inductance in such system can 
be characterized by the exponent vl = ln(125/33) / ln4 w 
0.96. On the other hand, the application of a simplified 
approach [18], which in terms of this work is equivalent 
to consideration of only (ALj )i, leads to a larger value 
ve = In5/ln4 w 1.16. That means that in the case of a 
Sierpinski gasket array the main contribution to the cor- 
rection to sheet inductance for small fields is coming from 
the term (AG n ) 2 , which can be associated with oscilla- 
tions of vortices. The results of this work show that in 



4 



periodic arrays the situation is qualitatively different and 
the main contribution to correction to inductance is com- 
ing from the suppression of the superfluid density related 
with a non-uniform distribution of the order parameter 
phases induced by the presence of vortices. 
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